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Abstract 

The aeroelastic stability of flap 
bending, lead-lag bending and torsion of a 
bearingless circulation control rotor blade In 
hover Is Investigated using a finite element 
formulation based on Hamilton's principle. The 
fl exbeam, the torque tube and the outboard 
blade are discretized Into beam elements, 
each with fifteen nodal degrees of freedom. 
Quasisteady strip theory Is used to evaluate 
the aerodynamic forces and the airfoil 
characteristics are represented either In the 
form of simple analytical expressions or In 
the form of data tables. The unsteady aero- 
dynamic effects are Introduced approximately 
through dynamic wake induced inflow modeling. 
The nonlinear equations of motion are solved 
for steady blade deflections using an itera- 
tive procedure. The flutter solutu n is 
calculated assuming blade motion to be a 
small perturbation about the steady solution, 
and the normal mode equations are used to 
reduce the number of equations. A correlation 
study of analytical results with the experi- 
mental data is attempted for selected 
bearingless bTade configurations with conven- 
tional airfoil characteristics. Then stability 
results are obtained for circulation control 
bearingless configurations consisting of a 
single flexbeam with a wrap-around type 
torque tube and the pitch links located at 
both the leading edge and the trailing edge 
of the torque tube. The stability is exa- 
mined at various thrust levels and collective 
pitch settings. 

NOMENCLATURE 

a. reference lift curve slope 

r (5.7/rad) 

c blade chord, m 


d pitch link location from torque 

tube center, m 

D blade section drag force, N 

e d aerodynamic center offset from 

elastic axis, positive aft, m 

e center of mass offset from elastic 

9 axis, positive forward, m 

I b blade mass moment of Inertia about 

0 flap axis, N-m-sec z 

L blade section lift force, N 


L ,L ,L aerodynamic force per unit length 
" in u,v,w, directions, N/m 

m mass per unit length of blade, 

N-ser./m 2 

m reference mass per unit length, 

0 N-sec 2 /m 2 

M 5 blade section aerodynamic moment 

* about midchord, N-m 

M aerodynamic moment per unit length 

* about elastic axis, N 

n number of elements 

s constant defining blowing distri- 

bution 

R rotor radius, m 

t time, sec 

u,v,w f elastic displacements in the x,y,z 

directions respectively , m 


c^ blade section drag coefficient 

c ft blade section lift coefficient 

c m blade section moment coefficient 

.5 about midchord 

C T rotor thrust coefficient — — — 7 

T p(nR) *R 


mV 


u blowing momentum coefficient 

Vfcpv z c 
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U T 

V 


V 1 


x.y.z, 


blade section normal velocity, 
m/sec 

blade section inplane velocity, 
m/sec 

blade section resultant air 
velocity, m/sec 


Induced Inflow, m/sec 


jet velocity (blowing), m/sec 
pitch link displacement, m 
undeformed blade coordinates, m 
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a blade section angle of attack, rad 

B p blade precone angle, rad 

s Lock number ^a p cR7l b 

o air density, N-sec 2 /m 4 


6T.5V variation of kinetic and strain 

energies respectively 

5W virtual work done due to aerodyna- 

mic loads 

e blade pretwist, rad 

X rotor induced inflow ratio, v^/ftR 


£,n,c deformed blade coordinates non- 

dimensional i zed wrt R 

a solidity ratio, blade area/disk 

area 

$ elastic twist about elastic axis, 

rad 

$ geometric apparent twist about 

deflected elastic axis, due to 
coordinate transformations, rad 


* dimensionless time, ftt 


<j v ,u) w ,u) fundamental coupled rotating lead- 

* lag, flap and torsion natural 
frequencies respectively 

ft rotor blade angular speed, rad/sec 

t ratio of modal damping to critical 

damping 

n t torque tube center offset from 

elastic axis, positive forward, m 

fl exbeam center offset from 
elastic axis, positive forward, m 

Introduction 

A circulation control rotor (CCR) utilizes 
circulation control (CC) aerodynamics for nfel.t 
rotor blade design. A CC airfoil is typically 
of quasi-el liptic profMe with 'ounded trailing 
edge, and a thin jet of air is 'lown from a 
spanwise slot (Fig. 1). The air jet remains 
attached over the curved profile because of 
Coanda effect "alance of centrifugal force 
and suction pressure). In a CCR, the thrust 
vector can be controlled by modulation of 
blowing as w"11 as geometric pitch. With a 
CCR, a high thrust is possible at reduced tip 
speeds and also the hub design can be 
simplified because of elimination of cyclic 
pitch. The application of CC technology is 
currently being evaluated In the design deve- 
lopment of a full-scale rotor. One concern Is 
the Influence of blowing on the dynamics of 
the rotor blade. 


With the availability of Improved 
materials, recent rotor design trends are 
leaning towards hingeless blade con- 
figurations. A bearingless rotor Is one such 
example where flap and lag hinges as well as 
pitch bearing are eliminated, and these are 
replaced by a root flexure consisting of 
flexbeam(s) and a torque tube (Fig. 2). The 
torsional ly soft flexbeam(s) extends from the 
hub to about 15-40% of blade radius where it 
is connected to the main blade. The pitch 
control to the blade is applied through the 
torsional ly stiff torque tube by rotating it 
with the pitch link which elastically twists 
the flexbeam(s). This resultr in a multiple- 
load-path structure because of the redundancy 
of 1 ad paths at the flexure. This causes 
however a more involved dynamic analysis. 

The objective of the present paper is to 
examine aeroelastic stability of flap bending, 
lead- lag bending and torsion of a circulation 
control bearingless blade in hover. 

A general review on aeroelastic stability 
of a rotating blade with conventional aerody- 
namics is given in References 1-2. Chopra and 
Johnson 3 formulated and analyzed the flap-lag- 
torsion aeroelastic stability of a CCR blade 
in hovering flight. Tnree degrees of motion 
were considered: rigid flap, lag and feather 
rotations about hinges at the blade root. The 
CC airfoil characteristics were represented in 
the form of simple analytical expressions. It 
was shown that the trailing edge blowing can 
have a major influence on blade aeroelastic 
stability. Recently Chopra 4 analyzed the 
aeroelastic stability of flap bending, lead-lag 
bending, and torsion of a CCR blade in hover 
using a finite element formulation. The CC 
airfoil characteristics in the form of data 
tables were used. Results were presented for 
several hingeless blade configurations. Again, 
it was shown that the blowing has an Important 
influence on blade dynamics which must be con- 
sidered in rotor design. 

Sivaneri and Chopra 5 applied a finite ele- 
ment formulation to analyze the dynamics of a 
bearingless rotor blade in hover with conven- 
tional aerodynamic characteristics. Each of 
the flexbeams and the torqje tube were modelled 
as Individual elastic beams. The displacement 
compatibility conditions at the clevis, between 
the inboard flexure beams and the outboard 
blade, were satisfied. Results were also calcu- 
lated using a simple equivalent-beam modelling 
wherein a bearingless blade Is represented as a 
single beam with equivalent properties. 

Comparison of the two sets of results showed that 
the equivalent-beam modelling can be quite 
erroneous for some cases. 

In the present paper, the above finite 
element formulation Is modified to study the 
aeroelastic stability of a bearingless CCR 
blade in hover. The multibeams of the flexure, 
and the outboard blade idealized as an elastic 
beam, are all discretized Into beam elements, 
each element with fifteen nodal degrees of 
freedom. There is a continuity of axial 
displacement u, flap bending w and w 1 , lead-lag 
bending v and v 1 and geometric twist $ between 
elements. Quasisteady strip theory Is used to 
evaluate aerodynamic forces. The airfoil 
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characteristics are represented In the form of 
data tables. The Influence of unsteady aero- 
dynamics Is Introduced by using dynamic Inflow 
modelling. The formulation is made quite 
general for nonuniform blades keeping In view 
the application to different types of 
bearingless configurations. First, a correla- 
tion study of analytical results with the 
experimental data is attempted for selected 
bearingless configurations with conventional 
airfoil characteristics. Then, stability 
resul ty are obtained for circulation control 
bearingless configurations. 


The aerodynamic coefficients C i% C d and 
Cjn * 5 are taken from data tables, and the 
numerical values for these coefficients are 
given at small steps: Aa of lfe° and a C y of 

1/200. The C is blowing momentum coefficient 
defined as u 



where mVj is the jet momentum, ty£pV 2 is the 
dynamic pressure, and c is the blade chord. 


Formulation 

The formulation details can be seen In 
References 4-5. The blade 1$ treated as an 
elastic beam and undergoes axial displacement 
u, lead-lag bending displacement v, flap 
bending displacement w and elastic twist 4 
about a deformed elastic axis. Fig. 3 shows 
the deformed as well as undeformed blade posi- 
tions. The rectangular coordinate system 
x,y,z is attached to the undeformed blade, 
wherein the x-axIs coincides with the elastic 
axis. A point P on the undeformed elastic 
axis undergoes displacements u, v, w in the 
x,y,z directions respectively and occupies the 
position P* on the deformed elastic axis. 

Then the blade section containing point P* 
undergoes a rotation e. about deformed 
elastic axis. 

\ * 0 ♦ ♦ (la) 


4 « 4 - / v*Vdx (lb) 

0 

where e Is pretwist, 4 Is the geometric twist 
with respect to the undeformed axis 
(compatible with u,v,w) and 4 is the elastic 
twist about the deformed elastic axis £. The 
formulation Is based on Hamilton's principle 

to 

n$U - 6 T - 6 W) dt 3 0 (2) 


where 6 U, 6 T and 6 W are respectively the 
variation of strain energy, the variation 
of kinetic energy and the virtual work done. 
These energy expressions are made independent 
of the time derivatives of virtual displace- 
ments, <5u, 6 v, 6 w and 64 and hence Eq. (2) can 
be written as 

6 U - 6T - 6 W 3 0 (3) 

The aerodynamic forces are obtained using 
quasisteady strip-theory approximation. 

Forces of non-clrculatory origin are also 
Included. The section lift, drag and moment 
about the mid-chord (per unit span) are 

L oV 2 c C t (a, C y ) 

0 • * PV*c C d (a, C v ) (4) 

M .5-T pv * c \ 5 (a> V 


The wake induced inflow Is assumed uniform 
along the length of the blade and the steady 
component is calculated from the momentum 
theory 

X 0 - K h (C T /2)^ ( 6 ) 

where Is an empirical factor and is assumed 
to be 1.15 and * t is the steady thrust coef- 
ficient. With blade vibratory motion, 
unsteady flow environments are created and 
which will naturally result in dynamic Induced 
Inflow condition. For hover, a simple dynamic 
inflow model is used 

rUx-tL ,7) 

4Xg 

where t is the time lag in arc radian and can 
be approximately taken as .85 / Aq. The A is a 
time induced inflow component, a perturbation 
about the steady component Xq. The aCj is the 
perturbation thrust component caused by blade 
motions. The blowing momentum coefficient 
C y Is not uniform along the length of the 
blade, and a general distribution is used for 
the formulation 

c v - (8) 

where s is a constant, i 3 r/R and C yT is blowing 
coefficient at blade tip. For simplicity of ana- 
lysis, it Is assumed that the blowing coefficient 
is constant within each element (based on mid- 
point of element). 


Finite Element Discretization 


The finite element formulation is based on 
energy principles (Hamilton). The flexbeam(s), 
the torque tube and the main outboaru blade ire 
all discretized Into a number of beam elements. 
Each element (Fig. 4) consists of fifteen degrees 
of freedom. There 1$ a continuity of u, v, v', 
w, w* and 4 between elements, and there are three 
Internal nodes, two for u and one for 4 . The 
distribution for deflections over an element are 
represented in terms of element degrees of freedom 
and shape functions; a second order polynomial for 
4 and cubic polynomials for u, v and w. 

Hamilton's principle In discretized form for 
n elements is expressed as 

n 

J («u 1 - 4T r dup ■ 0 (9) 
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where 6Uj, 6Tj and 6Wj are virtual energy 

contributions from 1 th element. During the 
assembly of element matrices, one has to use 
the displacement compatibility conditions at 
Interelement boundaries to form global matr1 r 
ces. A simple analytical model of bearingless 
blade Is shown in Fig. 5. The displacement 
compatability at the clevis is 



i i i 

w t s w f = w b 


♦t = ♦f 3 ♦b 

where subscripts t, f and b respectively 
represent torque tube, fl exbeam and b’ade. 

Fig. 6 shows the root end of the torque 
tube with the control actuator located at the 
leading edge. The blade pitch is changed 
through the up and down movement of pitch 
link. The torque tube typically is very stiff 
torsional 1y as compared to the flexbeam. 

Therefore any vertical movement of the pitch 
link results in nearly rigid body pitch for 
the torque tube and an elastic twist distribu- 
tion for the flexbeam. To obtain different 
blade pitch angles, one has to adjust the ver- 
tical position of pitch link, W p . The pitch 
link flexibility is represented by spring 
stiffness K p . Because of the pitch link 
flexibility, there will be an extra strain 
energy contribution for the last element of 
the torque tube, say V p , 

v p " i K p t w i + +i « * V - v> 2 HD 

This will modify the element stiffness matrix 
and the load vector. The relationship between 
blade pitch and pitch link displacement w p Is 
calculated iteratively. 

The assembly of n elements yields the 
equation of motion in terms of nodal displace- 
ments fq} as 

[M(q)]{q}+ [C(q)]|q} ♦ iK(q)]{q) - {Q} ( 12 > 

These are nonlinear equations in q. The next 
step is to a*>ply geometric boundary conditions. 
For the torque tube end, the axial and lead-lag 
displacements are freely permitted and the pitch 
link joint Is pin-ended. Therefore, there Is 
no geometric constraints for the torque tube. 

On the other hand, the flexbeam Is cantilevered 
at the hub and therefore the displacements u, 
v, v*, w, w' , and * are all zero at the root 
of the flexbeam. The boundary conditions are 
applied to the global equations (12) by canceling 
out vhe rows and columns corresponding to these 
constraint displacements. 


Solution Procedure 


The first step is to determine blade steady 
equilibrium position. For a known thrust level, 
the collective pitch is calculated. For this 
collective pitch, an approximate value of pitch 
link position w p is determined. Now with a known 
Wp as a boundary condition, the blade steady 
deflected position is calculated from the nonlinear 
equations (12), after dropping time dependent 
terms. 

The second step is to obtain the coupled 
natural vibration (rotating) characteristics of the 
blade about its steady equilibrium position. Then, 
the flutter solution is obtained by assuming vibra- 
tory motion to be small perturbation about 
equilibrium position. For this, the normal mode 
equations are obtained and an eigen analysis is 
made. The nature of the complex eigenvalues 
explains whether blade is stable or not. 


Results and Discussion 


Numerical results are calculated for 
bearingless blade configurations, consisting of 
single flexbeam with wrap-around type torque tube 
(as shown in Fig. 2). For calculations, the blade 
is discretized Into seven elements; three elements 
for main blade, two elements for flexbeam and two 
elements for torque tube. The normal mode stability 
solution is calculated using six coupled rotating 
modes. For analysis, the flexbeam and the torque 
tube are modelled as individual beams. At the root 
end, the flexbeam is rigidly fixed, whereas for the 
torque tt be there are no constraints on displacements 
except for a spring restraint in the pitch link 
direction. A limiting case of rigid pitch link will 
result in complete displacement constraint at the 
pitch link location. 


First, a correlation study of analytical 
results with the experimental data is attempted for 
selected bearingless blade configurations with con- 
ventional airfoil characteristics. The experimental 
stability data for a model rotor is taken from Ref. 

6. The model rotor characteristics are Lock number 
y = 5.9, solidity ratio o * 0.03, three-bladed, and 
zero precone. The airfoil characteristics used are: 

C t « 0.15 + 5.73a 

C d = 0.0079 + 1.79a 2 

C m * -°* 012 

The idealized nondimensional structural properties 
used for elements are given in Table 1. Three differ- 
ent bearingless rotor configurations are considered 
depending on the pitch link location. The pitch links 
are located at a distance of 0.085R In the radial 
direction from the rotation axis. Fig. 7 presents 
the lag mode stability results for case I, where the 
pitch link Is located at the leading edge of the 
torque tube. This positioning of pitch link will 
cause a positive pitch-flap coupling (flap up 
causing nose down pitch). In the figure, the 
damping of the fundamental lag mode In the form of 
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real part of the eigenvalue Is plotted for dif- 
ferent collective pitch angles. For this lag 
mode, a structural damping Idealized by an 
equivalent viscous damping ratio of It Is used. 
Based on the experience of stability correla- 
tions of other authors. It appears that the 
present analytical results are in good 
agreement with the experimental results. Fig. 

8 shows the lag mode stability results for case 
II, where the pitch link Is located at the 
trailing edge of the torque tube. This will 
cause a negative pitch-flap coupling. In Fig. 
9, the results are presented for case III, 
where one pitch link Is located at the leading 
edge and the other at the trailing edge of the 
torque tube. This will not cause any pitch- 
flap coupling, but It will raise the torsional 
frequency of the blade. Again, the agreement 
of analytical results with the experimental 
data appears quite satisfactory. It is also 
noted here that the inclusion of dynamic Inflow 
has only a slight Influence on lag mode 
stability. 

The subsequent results are obtalred for 
oearlngless CCR blade configurations. These 
results are calculated for a CCR blade with 
Lock number y * 7.2 : solidity ratio a * 0.13, 
four bladed, and zero precone. For simple 
airfoil characteristics, the following analy- 
tical expressions are used: 

C. = 0.3 + 6.7a + 16. 1C * 67 
l y 

C, - 0.026 - 0.3C * 67 
d u 

C ra - 0.06 + 1.34a - 0.644C U ' 67 

and these are gross representations of the data 
in Ref. 7. For table aerodynamics, the air- 
foil characteristics of a typical CC airfoil 
with single trailing edge slot are used. The 
slot height to chord ratio is taken as 0.002 
and the airfoil thickness to chord ratio is 
taken as 0.15. The chordwise offsets of the 
center of mass, the aerodynamic center, and 
the tension center from the elastic axis are 
considered to be zero, and the elastic axis Is 
assumed to be at mid-chord position. The non- 
dimensional structural blade properties for 
different elements are given in Table 2. For 
stability results, the Inherent structural 
damping Is assumed to be zero for all modes. 
Results are calculated for the CCR bearingless 
configurations, case III only, for which the 
pitch links are located at the leading edge 
and the trailing edge of the torque tube. 

Fig. 10 shows for trim solutions In hover, 
the blowing momentum coefficient at the blade 
tip plotted as a function of rotor thrust for 
several collective pitch angles. It Is 
assumed that the blowing coefficient varies 
Inverse to the square of the radial position 
(C u ■ C v y/e 2 ). Results are calculated using 
• simple expressions as well as airfoil tables. 
With airfoil tables, an iterative procedure 
based on Newton-Raphson method is used to 
calculate the trim solution. There is a con- 
siderable disparity between the two results. 
Tnls Is understandable since simple 
expressions are only gross representation of 


the airfoil characteristics below stall, 
whereas the table data covers the complete 
range of angle of attack and blowing coef- 
ficient. 

Fig. 11 shows the damping ratio c for 
three modes with less damping as a function of 
thrust level for a case with zero collective 
pitch (e 0 = 0). These less damped modes hap- 
pen to be fundamental lead-lag, fundamental 
torsion and second lead-lag modes. The other 
modes are more damped for this as well as for 
subsequent cases and hence are not plotted. 

These results are calculated using the airfoil 
tables. The negative value of damping repre- 
sents the instability condition of a mode. 

The lag mode is stable at low thrust levels 
and becomes unstable at high thrust levels. 

The torsion mode and the second lag mode are 
moderately unstable. All three modes can be 
easily stabilized with the inclusion of small 
amount of structural damping in these modes. 

Most rotor designs inherently have enough 
structural damping to stabilize these levels 
of instability. A similar type of stability 
characteristics were observed in Ref. 4 for 
hingeless CCR blades with zero collective 
pitch. In fig. 12, the stability results are 
shown for the same blade configurations using 
the simple analytical expressions presented 
previously for the airfoil characteristics. 
Though this approximation is quite gross as seen 
in the trim solution (Fig. 8), the stability 
results are quite reasonably predicted. The 
inclusion of dynamic inflow has a slight 
influence on lag mode stability, in fact, a 
destabilizing effect at high thrust levels. 

Fig. 13 shows the stability results for a 
negative collective pitch of -10° using the 
table aerodynamics. For this pitch setting, 
one needs a larger amount of blowing to acheive 
certain thrust level. The fundamental lag mode 
becomes unstable at low thrust levels, and 
becomes quite stable at high thrust levels. In 
fact, at low thrust levels one needs a larger 
amount of damping to stabilize this mode. The 
effect of negative pitch on torsion and second 
lag mode is comparatively less, it stabilizes 
the torsion mode somewhat. The effect of posi- 
tive collective setting on blade stability is 
shown in Fig. 14. For this pitch, one needs a 
small amount of blowing to achieve certain 
thrust level. Here, the fundamental lag mode 
gets very stabilized. The influence on the 
other two modes is again small. 


Conclusions 

The aeroelastlc stability of a bearingless 
circulation control rotor blade in hover Is 
examined using a finite element formulation. 
Airfoil characteristics are represented In the 
form of simple expressions as well as in the form 
of data tables. The flexbeam, the torque tube 
and the main blade are modeled individually as 
elastic beams. Numerical results are calculated 
for a circulation control bearingless rotor 
configuration consisting of a single flexbeam 
with a wrap-around type torque tube and the 
pitch links located at both the leading edge 
and the trailing edge of the torque tube. 
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Based on this effort the following conclu- 
sions can be drawn; 

(1) A reasonable correlation of analytical 
results with the experimental data 

for selected bearingless configurations 
with conventional airfoils has introduced 
confidence to the methodology and the 
algorithms. 

(2) For flutter solution, one needs to include 
at least four modes; these are fundamental 
flap, lag, and torsion modes and the 
second lag mode. 

(3) Stability results have shown that the 
trailing edge blowing plays an important 
role in the determination of the aero- 
elastic stability of a rotor blade. 

(4) The inclusion of dynamic inflow has a 
slight influence on lag mode stability. 

(5) The negative pitch setting destabilizes lag 
mode, in particular, at low thrust levels. 

(6) The positive pitch setting stabilizes lag 
mode. 

(7) For a fixed thrust level, it appears that 
going to higher blowing is more destabi lizing. 

(8) The fundamental torsion and the second lag 
modes are weakly unstable for some thrust 
levels. The expected levels of internal 
structural damping appear adequate to 
stabilize these modes. 
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Table 1. 

Structu 

ral propertl 

es of elements 

for rotation 

speed of 

1100 RPM. 


Element 

Length 

Flapwlse 

Chordwlse 

Torsion 

Mass 

Torsion Inertia 



l/R 

EI/VV 

e w 2r4 

GJ/m o n 2 R 4 

m/m 0 

K X 


1 

0.367 

0.0055 

0.1501 

0.0029 

1.0 

o.ro9i 

Blade 

2 

0.367 

0.0055 

0.1501 

0.0029 

1.0 

0.00091 

Blade 

3 

0.069 

0.1216 

0.1216 

0.2433 

39.6 

0.0105 

Blade 

4 

0.113 

0.00158 

0.0052 

0.00021 

0.299 

0.000029 

Flex beam 

5 

0.085 

2.099 

2.0S9 

9.150 

72.6 

0.0346 

Flex bean 

6 

0.0564 

4.257 

4.257 

1,815 

7.63 

0.0020 

Torque tube 

7 

0.0564 

4.257 

4.257 

1.815 

7.63 

0.0020 

Torque tube 


Chord/radius * .0465 
K rt 

Pitch link spring * 171.06 

m 0 n 2 R 


Pitch link offset from torque tube center ^ * .0409 

a Of 3 0 


Table 2. Structural properties of elements for the circulation control bearingless blade. 


Element 

Length 

t/R 

Flapwlse 

Ely/V 2 *" 

Chordwlse 

El z /m 0 nV 

Torsion 

Gj/m 0 nV 

Mass 

m/m 0 

Torsion inertia 

kVr 2 

m 


1 

0.2 

0.0186 

0.2303 

0.0297 

0.7067 

0.000739 

Blade 

2 

0.2 

0,0372 

0.3938 

0.0557 

1.0 

0.000832 

Blade 

3 

0.2 

0.0929 

0.7133 

0.0929 

1.624 

0.001068 

Blade 

/ 

0.2 

0.1858 

0,2303 

0,00297 

1.383 

0.000099 

FI exbeam 

5 

0.2 

0.5573 

0.6687 

0.00297 

1.556 

0.000279 

FI exbeam 

6 

0.13 

0.0817 

0.5201 

0.1560 

1.398 

0.001397 

Torque tube 

i ' 

1 

0.13 

0.1486 
1 

0.3901 

0.2823 

1.549 

0.001366 

Torque tube 


Chord/radius * .1034 

K * 

Pitch link spring — — * 10000 \ 

m 0 fTR 

Pitch link offset from torque tube center ^ ■ .033 
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Fig. 5 Analytical model of a bearingless blade. 
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Fig. 6 Torque tube and pitch link. 
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Fig. 7 Lag damping as a function of blade I 

collective pitch, case I (pitch link j 

at leading edge of torque tube). * 
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Fig. 8 Lag Hamping as a function of blade 

colTective pitch, case II (pitch link 
at trailing edge of torque tube) 
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Fig. 10 Blade blowing coefficient as 
function of rotor thrust and 
collective pitch. 
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Fig. 11 Damping ratio for zero collective pitch 
(zero structural dampings, u =2.5, u =2.3, u sl7.4) 
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Fig. 12 Damping ratio for zero collective pitch 
(zero structural dampings, w v =2.5, w w =2.3, u *17.4) 
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Fig. 13 Damping ratio for negative collective pitch Fig, 14 Damping ratio for positive collective pitch 

(zero structural dampings, w *2.5, w *2.3, u *17.0) (zero structural dampings, « *2.", u *2.3, u *17.4) 
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DYNAMIC STABILITY OF A BEARINGLESS CIRCULATION CONTROL ROTOR BLADE IN HOVER 

Inderjit Chopra 


Jack Nielsen, NASA Ames Research Center : I have a serious question with regard to the funda- 

mental aerodynamics of the problem. You have assumed you can calculate the flutter using steady 
aerodynamic results. Now we have a boundary layer on the Coanda plate and it isn’t clear at all 
that at the frequencies of flutter that the boundary layer isn't going to have very important 
unsteady effects; we really don’t know. But I don't think we can assume the quasisteady assump- 
tion offhand. Letting the amplitudes get small doesn't get around this problem if the unsteady 
effects are coupled into the flutter. So I would be interested to know since we don't have any 
data how you can be sure that you have really solved the real flutter problem for the CC air- 
foil. 

Chopra : Are you posing the question [to me]? I do recognize that the unsteady aerodynamics for 

CCR is very important and it's not there. This is something that has to be looked at in the 
future. I am very much interested [and would] like to work on that prodem if you give us the 
money. 

Bill Warmbrodt, NASA Ames Research Center : The results that you showed for the circulation 

controlled airfoil showed damping ratios that were an order of magnitude lower than the hinge- 
less configuration that you first showed the results for. You made mention that structural 
damping, had It been included, would have stabilized some of those modis. I think everybody In 
the audience here has a pretty good feel for what the influence of structural damping Is for a 
hingeless rotor configuration. Would you say that you saw that same legree of sensitivity for 
the bearingless configurations that you analyzed? 

Chopra : I think your first question is . . . let's look at it this way. You get lower damping 

here than you see in the hingeless [rotor], I think that [you] have to go back and look first 
at the perturbation aerodynamics. Keep In mind that you are keeping your elastic axis at the 
half chord. Your perturbation aerodynamics has two components— [one] due to blowing, [one] due 
to conventional. The conventional part of the aerodynamics still has the quarter chord as the 
aerodynamic center and that will be destabilizing particularly for the torsion mode. So you 
expect co see lower damping. [This] means the people who have to design the X-Wing/CCR Rotor, 
they do need to — what should I say — add more structural damping or some sort of damping to 
stabilize these modes. I hope you may be agreeable to this point. 

Jing Yen, Bell Helicopter : I understand that you used the normal modes approach to solve the 

problem. Would you like to tell us where and how you obtained the normal modes for the redun- 
dant or multiple load path hub? Also how much confidence [do] you have in [the modes]? 

Chopra : The confidence is 100<. This is something [that is] a routine classroom problem. 

What \3 do is we look at first to the steady deflected shape and at that stage we solv* the 
global equations. Solving the global [algebraic, eigenvalue] equation is very routine these 
days. It doesn't take any extra time. But when you are trying to solve the complex eigenvalue 
problem it is no good to use, say, a 100 by 100 equation, but it is good to use a 6 by 6 
equation. There are two reasons for that--one is the computer and the second is the physical. 
You don't want to look at a hundred eigenvalues, you want to look at Just five or six and see 
what is happening to (the major] mode shapes. Reducing to normal modes is the same [as] if you 
had got the mode shapes using the Myklestad method and reduce them to normal modes. I don't see 
any difference from the . . . if you were to [model] a beam [using], say, the Myklestad approach 
and 20 elements, you may be using only four or five modes. So the procedure is Just the same. 
Did I satisfy? 

Yen : This was a finite element model? 

Chopra : It's a finite element model, yes. To start with it's a finite element model. W* get 

the natural mode shape of [the blade] using finite elements. Redundancy doesn't come int-j the 
picture anywhere. This is only the way you are arranging the equations, 

Peretz Friedmann : I think I misunderstood his question because I had the same concern about 
the redundancy. I tnink that when you have a redundant structure which is what [Jing] Yen 
alluded to you don't know exactly what the boundary conditions are and you might get incorrect 
mode shapes if you are not careful about the boundary conditions. I am not sure what those 
boundary conditions are because you have the cuff and you have the redundant load path and you 
really don't know what exactly the boundary condition is at the root. I think that's what he 
meant. 

Chopra : I think that [is true with] any problem. If we don't know the boundary conditions we 

can make an error in any analysis. Some of the configurations I've seen, the pitcn links seem 
to be . . . you can easily [get] the moment there and the only really important displacement Is 
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the vertical displacement. t3u* /e did make a parametric study where we tried to constrain the 
pitch link. We tried to put various types of constraints on the pitch link. It is not very 
sensitive on stability. 

W ayne Johnson : Aren’t your normal modes, though, calculated after you find the deflected solu- 

tion using th3 full finite element. 

Chopra : That’j right. 

Johnson : So the normal modes are found after you have identified the boundary conditions. 

Chopra : That’s right, aft identifying the boundary conditions, that's right. 

Johnson : I think modes in this case is simply a way to reduce the dimension of the state vec- 

tor. I don't think it really dot.3 anything more than that. 
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